Cosmologies of extended massive gravity 
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We study the background cosmology of two extensions of dRGT massive gravity. The first is 
variable mass massive gravity, where the fixed graviton mass of dRGT is replaced by the expectation 
value of a scalar field. We ask whether self-inflation can be driven by the self-accelerated branch 
of this theory, and we find that, while such solutions can exist for a short period, they cannot be 
sustained for a cosmologically useful time. Furthermore, we demonstrate that there generally exist 
future curvature singularities of the "big brake" form in cosmological solutions to these theories. 
The second extension is the covariant coupling of galileons to massive gravity. We find that, as in 
pure dRGT gravity, flat FRW solutions do not exist. Open FRW solutions do exist - they consist of 
a branch of self-accelerating solutions that are identical to those of dRGT, and a new second branch 
of solutions which do not appear in dRGT. 



INTRODUCTION AND OUTLINE 

An interacting theory of a massive graviton, free of 
the Boulwarc-Dcscr mode [1], has recently been discov- 
ered [2, 3] (the dRGT theory, see [4] for a review), al- 
lowing for the possibility of addressing questions of in- 
terest in cosmology. Pure dRGT massive gravity admits 
self- accelerating solutions [5-10], in which the dc Sitter 
Hubble factor is of order the mass of the graviton. Since 
having a light graviton is technically natural [11, 12], such 
a solution is of great interest in the late-time universe to 
account for cosmic acceleration. 

A natural question is whether a similar phenomenon 
might drive inflation in the early universe. To use the 
self-accelerating solution of massive gravity for inflation 
(i.e. "self- inflation" ) , the graviton mass would have to be 
of order the Hubble scale during inflation. Yet, we know 
that the current graviton mass cannot be much larger 
than the Hubble scale today [13]. 

Thus, for self-inflation to be possible, the graviton 
mass must change in time. One idea of how to realize 
this is to promote the graviton mass to a scalar field, <£, 
which has its own dynamics and can roll [9, 17]. The 
expectation value (VEV) of $ then sets the mass of the 
graviton. We can imagine that at early times $ has a 
large VEV, so that the graviton is very massive, and the 
universe self-inflates with a large Hubble constant. Then, 
at late times, $ rolls to a smaller VEV, self-inflation ends 
and the graviton mass attains a small value consistent 
with present day measurements. 

In this letter, we will see that, in practice, such an 
inflation-like implementation of massive gravity is diffi- 
cult to achieve in this model. Pure dRGT theory has 
a constraint, stemming from the Bianchi identity, which 
forbids standard FRW evolution in the flat slicing [9] (the 
self- accelerating solutions are found in other slicings). 
There appears an analogous constraint in the variable 
mass theory, and this constraint, while it no longer for- 
bids flat FRW solutions, implies that self-inflation cannot 



be sustained for a cosmologically relevant length of time. 
In addition, we show that non-inflationary cosmological 
solutions to this theory may exhibit future curvature sin- 
gularities of the "big brake" type. 

In the second half of this letter (which can be read 
independently from the first), we consider the covariant 
galilcon extension of massive gravity introduced in [14]. 
This theory has a new scalar degree of freedom tt, which 
describes brane bending in an additional spatial dimen- 
sion. (Unlike the variable mass theory, ir here does not 
act to set the graviton mass, which is fixed. So we are 
not interested in self-inflation here, but are just exploring 
the basic cosmological equations.) 

We derive the background cosmological equations for 
this theory, and find that the presence of the scalar leads 
to a more complicated constraint than in pure dRGT. We 
discuss the possible solutions in the case of zero and neg- 
ative spatial curvature. We find that, as in pure dRGT 
theory, this constraint forbids flat FRW solutions. For 
an open FRW ansatz, however, solutions can exist and 
they come in two branches. The first branch consists of 
self- accelerating solutions that are identical to the self- 
accelerating solutions of pure dRGT theory. The second 
branch consists of novel solutions which are not found in 
pure massive gravity. 



VARIABLE MASS MASSIVE GRAVITY 



We start with variable mass massive gravity. This is 
dRGT theory in which the graviton mass squared is pro- 
moted to a scalar field $, 



5 = 5, 



EH 
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where 



Seh — 



-M 2 / d^xJ—gR, 



(2) 



Inserting (6) and (9) into the action, we obtain the 
mini-supcrspacc action 



M 2 J d 4 x $ (^2 + a 3 C 3 + Q-iU) , (3) S R = 3Af£ / dt 



- / dV 



,g($)(a$) 2 + y($) 



(4) 



Here 0:3, are the two free parameters of dRGT theory. 
We have allowed for an arbitrary kinetic function g(&) 
and potential V(Q), so that there is no loss of generality 
in the scalar sector. The mass term consists of the ghost- 
free combinations [3], 



C 2 
C 3 



i([/C] 3 -3[/C][/C 2 ] + 2[/C 3 ]) 
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([/C] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]) 



where K% = 5% — \fg^ 3 r\^, 1 g^v is the non-dynamical 
fiducial metric which we have taken to be Minkowski, 
and the square brackets are traces. To work in the gauge 
invariant formalism, we introduce four Stiickelberg fields 
(j) a through the replacement g^ — > d^d^gab- 

Variable mass massive gravity was first considered in 
[9], and further studied in [16-19] (see also [15] for a more 
symmetric scalar extension of dRGT) . dRGT gravity has 
been demonstrated to be ghost-free through a variety of 
different approaches [20-24] , and the introduction of the 
scalar field does not introduce any new Boulware-Dcser 
like ghost degrees of freedom into the system [16]. 

For cosmological applications we take a Friedmann, 
Robertson- Walker (FRW) ansatz for the metric, so that 



ds 2 



-N 2 {t)dt 2 + a 2 [t)ni jdx l dx\ 



where 



1 — nr 2 



(0) 



(7) 



is the line element for a maximally symmetric 3-space of 
curvature k and r 2 = x 2 + y 2 + z 2 . We also take the 
assumptions of homogeneity and isotropy for the scalar 
field, 



$ = $(t). 



(8) 



Consider first the case of fiat Euclidean sections (k = 
0) . We work in the gauge invariant formulation, and the 
Stueckelberg degrees of freedom take the ansatz [9, 10]. 



4>° = /(*), 



(9) 



where f(t), like a(t), is a monotonically increasing func- 
tion of t. 



S mass =3M 2 I di$ 
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NF(a) - fG(a) 



dt a 3 
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(10) 

(11) 
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where 



F(a) = a(a - l)(2et - 1) 

+ ^(a-l) 2 (4a-l) + ^(a-l) 3 , 

G(a) = a 2 (a - 1) + a 3 a(a - l) 2 + ^±(a - 1) 



(13) 
(14) 



This mini-supcrspace action is invariant under time 
reparametrizations, under which / transforms like a 
scalar. 

There are four equations of motion, obtained by vary- 
ing with respect to F, N, $ and a. As in GR, the Nocthcr 
identity for time reparametrization invariance tells us 
that the acceleration equation obtained by varying with 
respect to a is a consequence of the other equations, so 
we may ignore it. After deriving the equations, we will 
fix the gauge TV = 1 (this cannot be done directly in the 
action without losing equations). 

Varying with respect to / we obtain the constraint 
pointed out in [17], 



C 



G(a) 



(15) 



where C is an arbitrary integration constant. (Note that 
the analogous equation in the fixed mass theory implies 
that a = constant, so there are no evolving flat FRW 
solutions in that case [9].) Varying with respect to N, 
we obtain the Friedmann equation, 



3M4 



H 



$F(a) 



= \g(m 2 + v, (16) 



and varying with respect to $ we obtain the scalar field 
equation 



$ + 3iT$ 



\g'(<S>)& + V ($) 
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Rather than solving the coupled second-order Einstein- 
scalar equations of motion, one can instead reduce the 
system to a single first-order Friedmann equation. The 
relation (15) can be used to eliminate $ and its first 
derivative from (16), which then becomes a first-order 



3 



differential equation in a which determines the scale fac- 
tor, 



H 2 = 



Q )lf2 \ni„ ( C \ G'(a) 2 a 2 
6M p - g I qjjj I G(a) 4 



(18) 



Once we have solved for the scale factor, the scalar $ 
is determined from (15) and the Stueckelberg field / is 
determined by solving (17) . 



Singularities 

One feature of this model that has not been noticed 
previously is that it allows for the possibility of curvature 
singularities at finite values of a. These can happen when 
the evolution attempts to pass through values of a for 
which the denominator of the right hand side of (18) 
goes to zero. 

When this happens a is finite, but the Hubble param- 
eter, and hence a, is blowing up. The scalar curvature 
is also blowing up, so this is a genuine curvature singu- 
larity; a "big brake" where the universe comes to some 
finite scale factor and gets stuck [25, 26]. Similar types 
of singularities also occur in DGP [27]. 

For example, Taylor expanding the denominator of 
(18) for large a we obtain a critical value of the scale 
factor at which the Hubble parameter diverges, 



= V3 



9(0) 
2A/2 



C 



3 + 3(^3 + «4 



1/3 



(19) 



Self-Inflation 

We now consider the possibility that the graviton has 
a large mass in the early universe, through some natu- 
ral displacement (and resulting VEV) of the scalar field 
from its true minimum near zero. We seek dynamics 
such that the scalar field slowly rolls down its potential, 
during which time the graviton remains massive, result- 
ing in a large self-acceleration of the universe. This self- 
acceleration comes from the second term on the left-hand 
side of the Friedmann equation (16). For this to be true 
self- acceleration this term should be much larger than 
both the scalar kinetic energy and potential on the right 
hand side, so that the acceleration is primarily driven 
by the modification to gravity and not by the scalar 
field. After many e-folds, $ should roll towards zero, 



1 Note that in general the Stueckelberg field cannot be chosen ar- 
bitrarily as in [17] but is non-trivially constrained by the choice 
of mass term, or in this case, kinetic function g(&). 



self-inflation should end, and the graviton mass should 
become small at late times. 

Thus, assume we have an inflationary solution a ~ 
e Ht , with H ~ constant. The scale factor is growing 
exponentially, so we Taylor expand the entire right hand 
side of (18) for large a, as 



V(0) 
3M P 



v'(o) 



(6 + 4a 3 + a 4 ) 



3 + 3tt3 + «4 



(20) 

We see that the dependence on all of the massive grav- 
ity modifications redshifts away exponentially, at least as 
fast as a -3 , and we are left with inflation driven only by 
the value of the potential at $ = 0. (In particular, con- 
tributions sensitive to the function g only start to enter 
at 0(l/a 7 ).) 

Said another way, we only have self-inflation if the 
quantity in (16) is approximately constant when 

a ~ e Ht . But the constraint equation (15) makes this im- 
possible: we see from (15) that $ ~ -i^, since G(a) ~ a 3 
for large a. Since F(a) ~ a 3 , the quantity ^^"^ behaves 
like ~ $ ~ a -3 , so it goes to zero exponentially fast and 
we cannot sustain self-inflation. 

Having encountered an obstacle to the possibility of 
self-inflation in the flat slicing, we now investigate the 
possibility in the open slicing (k < 0). Following [10] we 
take the Stueckelberg ansatz to be 



The mini-superspace action then becomes 



Sr 


= 3A4 J 


mass 


= 3A4 J 




= J dta 3 



N 



nNa 



NF(a,f)-fG(a,f) 



-7V~ 1 .g($)$ 2 - NV{$) 



(21) 

(22) 
(23) 
(24) 



where 



F(a, f) = a(a - (2a - v 7 ^/) 




G(aJ) = a 2 (a- 

<">4 




R/) 2 (4a- 



nf) + a 3 a(a 



a 4 

T 

?\2 



(a 




(26) 



Again, we have time reparametrization invariance so we 
can ignore the acceleration equation, and we will fix the 
gauge N = 1 after deriving the equations of motion. The 
constraint equation arising from varying with respect to 
/is 



(a- 



dG(a,f) 
da 



+ G(a,/)$ = 0. (27) 
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The Fricdmann equation obtained by varying with re- 
spect to N is 



H 



\g(<S>)tf + (28) 



and the equation of motion for $ is 



$ + 3i?3> 



i 5 '($)$ 2 + y'($) 



3A« 



F(q,/) ; G(o,/) 



(29) 



In order to obtain inflation driven by the graviton 
mass, the term $>F(a, f)/a 3 in the Friedmann equation 
(28) must be approximately constant when a ~ e Ht . Re- 
arranging the constraint equation (27) to isolate $ gives 



4> 
4> 



H 



da 

a J G(a,f) 



(30) 



Now there are three possibilities: the first is that 
a(t) ~ e Ht grows faster than f(t). In this case, the right- 
hand side of (30) approaches a constant at late times, 
<j>/$ — 3i/ + 0(1/ a), which tells us that $ decreases 



exponentially at late times, ~ e 



This, in turn, 



implies that the self-acceleration quantity $F(a, f)/a 3 
in (28) decreases exponentially like since F(a,f)/a 3 
approaches a constant. So once again, we cannot sus- 
tain inflation in this case. The second possibility is that 



f(t) grows faster than a(t) 



In this case, the right- 



hand side of (30) goes to zero at late times, so $ becomes 
constant. The self-acceleration quantity $F(a, f)/a 3 in 
(28) now grows without bound as ~ / 3 (— «) 3 / 2 /a 3 at 
late times, so again we do not achieve sustained infla- 
tion. Finally, there is the possibility that f(t) ~ e~ Ht , 
growing at the same rate as a(t). This case follows the 
same pattern as the first possibility - the right-hand side 
of (30) approaches a constant at late times, and the self- 
acceleration quantity decreases exponentially. 

In summary, flat and open FRW solutions exist in 
mass-varying massive gravity, but the constraint equa- 
tion (the one which forbids flat FRW solutions in 
pure dRGT theory) does not allow for long-lasting self- 
inflation. Finally, homogeneous and isotropic closed 
FRW solutions are not possible for the same reason they 
are not in dRGT - the fiducial Minkowski metric cannot 
be foliated by closed slices. 



THE DBI MODEL 

We now turn to the cosmology of the coupled galileon- 
massive gravity model introduced in [14]. One way to 
construct the theory is to consider a dynamical metric, 
g^vi on the world volume of a 3-brane in a five dimensional 



flat bulk, with embedding coordinates X A (x), coupled 
to the induced metric through dRGT interactions. The 
worldvolume theory is 



S[g, g] = Sn[g] + S mass [g,g] + Sir[g] 



(31) 



where Sr and Smass are the usual Einstcin-Hilbert and 
dRGT mass terms except that now 



g^ = d^X A d„X B VAB 



(32) 



is the pull-back of the 5D Minkowski metric to the 4D 
brane. The action S„ consists of ghost-free scalars con- 
structed from g^y and the extrinsic curvature [28]. The 
only term which will be relevant for us is the DBI term 



-A 4 



xv -g, 



(33) 



where A is some characteristic mass scale. The other 
possible terms in all vanish on the cosmological back- 
grounds we consider, so we may exclude them without 
loss of generality. The embedding functions are scalar 
fields; the first four become the Stueckelberg fields, and 
the fifth becomes an extra physical scalar, the brane 
bending mode 7T, 



X' 



(34) 



Consider first the case of flat Euclidean space (k = 0). 
We take the FRW ansatz (6) and (9), along with ho- 
mogeneity for the scalar tt = n(t). The mini-superspace 
action is 



S R = 3AI£ I dt 



"Smass — 3AI 

S„ = -A 4 



p I dtm z 



N 



NF(a)-G(a)^f 2 -ir 2 



^2 



(35) 
(36) 
(37) 



where F(a) and G(a) are as in (13) and (14). 

The mini-superspace action is invariant under time 
reparametrizations, under which both / and 7r transform 
like scalars. Thus, as before we may ignore the accelera- 
tion equation of motion obtained by varying with respect 
to a, since it is redundant. We fix the gauge N = 1 after 
deriving the equations. 

Varying with respect to N we obtain the Fricdmann 
equation 



0. 

The / and tt equations are, respectively 

/ 



SS _ d 
Jf~di 



(ZMlm 2 G(a) + A 4 



P 



(38) 



0, (39) 



5 



6S_ 

Stt 



d 

"dl 



(3Mpm 2 G(a) + A 4 ) 



P - * 2 



0. (40) 



Taking the combination /fj + ^ff of the two equa- 
tions of motion (39) and (40) , we arrive at the following 
constraint equation, 



P-n*-[G(a)}=0. 



(41) 



The square root cannot be zero because it appears 
in the denominator of (39) and (40), so we must have 
^ [G(a)] = 0, which implies that a must be constant and 
thus there are no evolving flat FRW solutions. This is the 
same phenomenon as in pure dRGT theory [9], and we 
see that the DBI extended theory has a similar structure. 
The inclusion of matter minimally coupled to does 
not affect this conclusion because it does not contribute 
to the equations (39), (40) and hence the constraint (41) 
is unchanged. 

We now turn to the case of nonzero spatial curvature 
k < 0. Using again the ansatz (21) we obtain 



Sr = 3Mp / di 



N 



nNa 



(42) 



Smass = 3M P / di rrv 



NF(a,f)-G(a,f)Jp-^ 



Stt = —A / dtV/ 2 -^ 2 , 



(43) 
(44) 



with F(a,f) and G(a, f) as in (25) and (26). 

Varying with respect to N and setting N = 1, we ob- 
tain the Fricdmann equation 



H z + — + m 



F(aJ) 



= 0. 



(45) 



The / and tt equations are, respectively 



SS_ 
Sf 



3M|m 



8F(aJ) dG(aJ) 



d 
di 



d.f df 
(3M^m 2 G(a, f) + A 4 ) 



f 2 -n 2 



0. 

(46) 

= 0. 
(47) 

Taking the combination /|j + 7i"ff of the two equa- 
tions of motion (46) and (47), and using the relation 



5S_ 

Stt 



d 
df 



(3M 2 m 2 G(a,f)+A 4 ) 



r- 



dF(a.f) _ 
df 



: 9G(aJ) 
" da i 



which can be checked straight- 
forwardly from the definitions (25) and (26), we arrive at 
the following constraint equation, 



dG(aJ) 
da 



0. 



(48) 



There are now two possible branches of solutions. The 
first consists of the solutions for which ^p- = 0. Solv- 
ing this algebraically for /, we find / = C± (-^J for some 
constant C± related to the coefficients ct2 , (there can 
be two branches here, since we must solve a quadratic 
equation for /). Then, reinserting this into (45), we ob- 
tain a modified Fricdmann equation. The modification 
F ^ a 'a? 1S a cons tant, depending only on a 2 , a^, so that 
when this constant is negative, we have self-acceleration 
with /9 so if ~ Mpm 2 . These solutions are exactly the same 
self- accelerated solutions that exist for the pure dRGT 
theory [10]. The solution for tt can then be determined 
by solving (47). 

However, for the theory at hand, there exists a new 
possibility. This second branch consists of solutions for 
which 



aV/ 2 



(49) 



In the case of the pure dRGT theory where tt = 0, this 
branch gives only solutions for which a = ^—nt, which is 
just Minkowski space in Milne coordinates. Here we have 
the possibility of non-trivial solutions on this branch. 
Solving (49) for tt gives tt = ±/^/l + p-, and substitut- 
ing this into the tt equation of motion (47) we see that 
/ cancels and we are left with an algebraic equation in / 
which can be solved, 



A 4 



/ 



3w 2 M|a 3 



C 



3m 2 M 2 a 2 y 



(50) 



where C is the integration constant from integrating (47). 
Using this to eliminate / from the Friedmann equation 
(45) yields an equation of motion for the scale factor, one 
which is a fourth order polynomial in H , and which can 
have solutions with non-trivial evolution. 



SUMMARY AND CONCLUSIONS 

In this letter we have explored the cosmology of two 
different extensions of the dRGT massive gravity theory. 

The first model is variable mass massive gravity, which 
results from promoting the fixed mass term of the dRGT 
model to the vacuum expectation value of a dynamical 
scalar field, as suggested in [9]. In dRGT theory, there is 
a constraint equation that forbids non-trivial flat FRW 
solutions, though self-accelerating open solutions exist. 
In the variable mass theory, the form of the constraint 



6 



is different, and it no longer forbids fiat FRW solutions. 
The constraint, however, makes it difficult to realize the 
idea of self-inflation, i.e. using the self-acceleration prop- 
erties of massive gravity in the early universe to drive 
inflation. Furthermore, we have demonstrated for the 
first time that a large class of cosmologies within these 
models exhibit a future curvature singularity of the "big 
brake" form. 

The second extension is the coupled galileon-massive 
gravity scalar-tensor theory of [14]. In the flat slicing, 
the constraint equation takes a similar form to that of 
the pure dRGT theory, and it similarly rules out non- 
trivial FRW solutions. In the open slicing, the theory 
reproduces the self-accelerating branch discovered in pure 
massive gravity, and in addition, provides a new branch 
of evolving cosmological solutions, the detailed general 
properties of which will be the topic of future work. 
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